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Abstract 

We have studied, with the KLOE detector at the DA$NE ^-Factory, the dynamics 
of the decay rj — ^ tt"^ tt" tt^ using data from the radiative $ — ^ 777 decay for an 
integrated luminosity L = 450 pb~^. Prom a fit to the Dahtz plot density distribution 
we obtain a precise measurement of the slope parameters. This should allow to 
improve the knowledge of the decay amplitude which is sensitive to the u-d quark 
mass difference. We also present new best results on the C-violating asymmetries 
in the 77 — tt"*" 7r~ vr'^ decay. 
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1 Introduction 



The decay of the isoscalar r] into three pions occurs through isospin violation and thus 
is sensitive to the up-down quark mass difference. The electromagnetic corrections 
to the decay are small (Sutherland's theorem [1]) and do not affect significantly the 
rate or the Dalitz plot density [2]. 

Neglecting electromagnetic corrections, the decay amplitude is given by [3]: 



where 



2 "2 

mt — m 



Q' = ^1—^ (2) 

is a combination of quark masses and m = | (m„ + rrid) is the average u, d quark 
mass, = 92.4 MeV is the pion decay constant and M{s, t, u) contains all the 
dependence of the amplitude on the Mandelstam invariants. Since the decay rate is 
proportional to Q^"^, 

r (r/ ^ 7r+7r-7r° ) oc \ A\^ oc (3) 
the transition — 37r is, in principle, a very sensitive probe to determine Q. A 
theoretical estimate of Q is obtained from the Dashen theorenill [4] according to 
which 



^Dashen ™2™2 _ 2 , 2 _ 2 l^^'J^J' 

rrij^o ^ ''''n+ ''^Tr" 

Of course, in order to extract the quark mass ratio from the decay width, one needs 

an accurate description of M{s,t,u). 

At lowest order 

3s — Am ^ 

M{s,t,u) = '-^^^. (5) 
— 

a well known result, based on Current Algebra. Integrating eq.(5) over the phase 
space and using the quark masses as estimated by Leutwyler [5] one obtains the 
following prediction for the decay rate [3] : 

ytheo (^^ _^ 7r+7r-7r° ) = 66 eV (6) 
which strongly disagrees with the experimental value [6]: 

T^^'P [t] 7r+7r-7r° ) = 295 ± 16 eV. (7) 
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The Dashen theorem states that in the chiral limit the charged kaon and pion electro- 



magnetic mass shifts are the same: — m^o = [^"k^ ~ n^\o 
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A one-loop calculation within conventional chiral perturbation theory (xPT) [7], 
improves considerably the prediction: 

r'^^" (t] 7r+7r-7r° ) ~ 167 ± 50 eV. (8) 

but is still far from the experimental value. 

Higher order corrections discussed in the literature [3] contribute to increase the 
theoretical prediction, but cannot fully account for the discrepancy with the data. 
More recently a good agreement has been found combining the effective chiral La- 
grangian with a non perturbative scheme based on coupled channels and Bethe 
Salpeter equation [8]. 

A significant violation of the Dashen theorem could in principle account for the dis- 
crepancy; however it should be claimed only after demonstrating agreement between 
theory and experiment on the M{s, t, u) amplitude behaviour over phase space. 
The above discussion therefore motivates a precise measurement of the ^ 37r dy- 
namics through the study of the Dalitz plot. 

The Dalitz plot of the t] t^^-k^-k^ decay is described by two kinematic variables. 
By convention these are defined in terms of the kinetic energies of the pions T^_, T_ 
and To in the r] rest frame: 

Qr) — rrirf — 2m^+ — 171^^0 , 



where X and Y respectively vary in the range [—1, 1] and [—1, 0.895]. 

The decay amplitude is then expanded around the center of the Dalitz plot (X=Y=0) 

in powers of X and Y as: 

\A{X, F)|2 ~ 1 + ar + bY^ + cX + dX^ + eXY + .... (9) 

and the parameters {a,b,c,d,e, ... ) of the expansion are fitted to the experimen- 
tal data. Any odd power of X contributing to |74(X, would imply violation of 
Charge Conjugation. 
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2 The KLOE detector 



Data were collected with the KLOE detector at DA$NE [9], the Frascati e^e~ 
collider, which operates at a center of mass energy ^/s = M$ ~1020 MeV/c^. The 
electron and positron beams collide with a crossing angle of tt — 25 mrad, resulting 
in a small momentum ( ~ 13 MeV/c in the horizontal plane) of the produced 
mesons. 

The KLOE detector consists of a large cylindrical drift chamber (DC), surrounded by 
a fine sampling lead-scintillating fibers electromagnetic calorimeter (EMC) inserted 
in a 0.52 T magnetic field.. 

The DC [10], 4 m diameter and 3.3 m long, has full stereo geometry and operates 
with a gas mixture of 90% Helium and 10% Isobutane. Momentum resolution is 
< 0.4%. Position resolution in r — is 150 /xm and cx^ ~ 2 mm. Charged tracks 
vertices are reconstructed with an accuracy of ~3 mm. 

The EMC [11] is divided into a barrel and two endcaps, for a total of 88 modules, 
and covers 98% of the solid angle. Arrival times of particles and space positions of 
the energy deposits are obtained from the signals collected at the two ends of the 
calorimeter modules, with a granularity of ~(4.4 x 4.4) cm^, for a total of 2240 cells 
arranged in five layers. Cells close in time and space are grouped into a calorimeter 
cluster. The cluster energy E is the sum of the cell energies, while the cluster time 
T and its position R are energy weighted averages. The respective resolutions are 
aE/E = 5.7%/^E{GeV) and ar = 57 ps/^E{GeV) © 100 ps. 

The KLOE trigger [12] is based on the coincidence of at least two energy deposits 
in the EMC, above a threshold that ranges between 50 and 150 MeV. In order to 
reduce the trigger rate due to cosmic rays crossing the detector, events with a large 
energy release in the outermost calorimeter planes are vetoed. 



3 Signal selection and efficiency 

At KLOE T] mesons are produced in the process — > 777 , so to study the dynamics 
of 77 — > 7r''"7r~7r° the final state 7r+7r^777 is used to which corresponds a BR value 
of: 

BRtot = BR{(p^ri-f)xBR(ri^ 71^%-%° ) x BR (tt^ ^ 77) = 2.9 ■ 10"^ 
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There is no combinatorial problem in the photon pairing because this decay chain is 
characterised by an energetic monochromatic recoil photon, with E^^^^ ~ 363 MeV, 
quite well separated from the softer photons from tt" decay. 

In the data analysis a photon is defined as a cluster in the EmC with no associated 
track detected in the Drift Chamber (DC) and with |(T — ^)| < 50"^; where T is 
the arrival time on the EmC, R is the distance of the cluster from the vertex, c is 
the light speed and ar is the time resolution. 
The events selection is performed through the following steps: 

(1) Events are selected starting from a very loose offline preselection consisting of a 
machine background filter (FILFO) and an event selection procedure (EVCL) 
that assigns events into categories [13]. 

(2) One charged vertex is required inside the cylindrical region r < 4 cm, |2;| < 8 cm 
and 3 photons with 21° < 9^ < 159° and > 10 MeV. The probability of a 
photon to fragment in more than a cluster (splitting) is reduced by requiring 
the opening angle between any pair of photons to be > 18°. 

(3) E ^7 < 800 MeV. 

(4) A constrained kinematic fit is performed imposing 4-momentum conservation 
and t — — for each photon, without imposing mass constraint both on rj and 
tt". 

A cut on the probability is done, P(x^) > 1%- 

This fit improves significantly the resolution on photon energies. The distri- 
bution of the fit shows a satisfactory Data-MC agreement as shown in fig.l: 




Fig. 1. distribution for the kinematic fit. Left: linear scale; Right: log scale. 

(5) Finally we require: 

• 320 MeV < E^rec < 400 MeV for the recoil photon (to reduce residual back- 
ground from 4> —>■ Ks Kl events) . 
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• E^++E^- < 550 MeV (to reduce residual background from — * 7r"'"7r~7r°events). 

• the invariant mass of the two softest photons: M^^ G [110, 160] MeV (to re- 
duce residual background from rj — > 7r"'"7r~7r° decays with 7r° e'^e~j and 
from (f) ^ T]^ events with rj 7r"'"7r~7). 

The selection efficiency is determined by using the KLOE MonteCarlo (MC) simu- 
lation program [13] and checked on data by means of control samples. In particular: 

• The trigger efficiency evaluated by MonteCarlo is 99.9%, and we have excellent 
Data-MC agreement for the trigger sectors multiplicities. 

• The effect of the event classification procedure (EVCL) and machine background 
filter is evaluated using a downscaled sample of non filtered data to which we 
apply much less stringent cuts in order to have a "minimum bias" selection. On 
signal events the efficiency of the minimum bias selection is 99.88%. 

We have found that the EVCL procedure introduces a signal loss of ~ 1.5%, ex- 
pected also from Monte Carlo estimates. The corresponding bias on the Dalitz 
plot parameters measurement has been evaluated and included in the systematic 
error. No bias is introduced by the FILFO procedure. 

• The tracking and vertex efficiencies have been estimated from the Data-MC ratio 
observed for a control sample of — > 7r+7r~7r° events selected to have charged-pion 
momenta in the same range as those from the rj 7r"'"7r~7r° decay [14]: 



The Data-MC ratio of efficiencies is flat all over the momentum spectrum, thus 
introducing no bias in the Dalitz plot fit. We recall that all the variables used 
in the fit are evaluated in the rj rest frame, which is boosted with respect to the 
laboratory by about 363 MeV/c momentum. Therefore to each momentum bin in 
the lab frame corresponds a much wider momentum interval in the 77 cms: any 
Data-MC discrepancy in the lab frame is further diluted by this effect. 

• A correction to the detection efficiency for low energy photons has been obtained 
by comparing the photon energy spectrum of a data subsample to the expected 
MC spectrum; the average correction factor is 0.964. 

The overall selection efficiency, taking into account all the Data-MC corrections is 



i^TRK^VTx) 



Data 



MC 



= 0.974 ± 0.006. 



(10) 
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found to be e = (33.4 ± 0.2)%. The expected background contamination, obtained 
from MC simulation is as low as 0.3%. 



After the background subtraction the observed number of events is: 

Nobs = 1.337 ± 0.001 Mevts. 
The distribution of the data in the Dalitz plot is shown in fig. 2. 




Fig. 2. Dalitz-plot distribution for the whole data sample. The plot contains 1.34 millions 
of events in 256 bins. 

The signal selection efficiency e {X, Y) as function of Dalitz plot point is obtained 
by MC, for each {X, Y) bin, as the ratio: 

where Nrec {N, Y) and Ng^n {N, Y) are respectively the reconstructed and generated 
Dalitz distributions. This definition of efficiency takes into account the smearing 
effects due to the finite resolution in X, Y. This approach is equivalent to the use of 
the complete four- dimensional smearing matrix as long as the MC correctly repro- 
duces the Dalitz plot shape. To this aim we have used a data subsample to obtain 
a first estimate of the Dalitz plot parameters which we have used in final MC. 
The efficiency e {X, Y) has a smooth behaviour all over the Dalitz plot. The pro- 
jections of e {X, Y) on the X and Y axis are plotted in fig. 3. While the efficiency 
appears to be rather fiat on X (and symmetric in X as expected), it decreases in an 
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Fig. 3. Left: Efficiency versus X. Right: Efficiency versus Y. 



approximately linear way with Y. In fact a large value for Y means a low- momentum 
TT^ in the decay to which corresponds a lower tracking/vertexing efficiency. 
The expected resolutions on the Dalitz variables {X, Y) are shown in fig. 4. The core 
resolution on X is about 0.02, due to our excellent momentum resolution for charged 
tracks. The Y variable, which is proportional to the kinetic energy, is evaluated 
[15] as the average between the "direct" determination obtained from the energy and 
direction of the two clusters associated to the 7r° — > 77 decay and the "indirect" 
determination : Tq — — {Et^+ + Et^-) — M^o . This leads to a core resolution on 
Y of about 0.02. 
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Fig. 4. Resolutions on X (left) and Y (right) according to MC. The curves are fitted to a 
sum of four gaussians. 



9 



4 Fit of Dalitz plot 



The fit to the Dahtz plot is done using a least squares approach. Let |^(^, be 
the theoretical squared amplitude: 

\A{X, Y)\^ = A^(l + aY + bY^ + cX + dX^ + eXY + ...). (12) 

with N being a normalization constant. 
Then the is defined as: 

^ ? ? I ^, j '''' 

where dPh(X,Y) is the phase space element and for each bin 

• Nij is the number of observed events in the bin, 

• Eij is the signal selection efficiency, 

• {xi, Xi + Ax) and {i/j, i/j + Ay) are the bin boundaries. 



All the bins are included in the fit except those intersected by the Dahtz plot contour. 
The fit procedure has been tested on Monte Carlo and has been verified to correctly 
reproduce in output the input values of the parameters. 

When applying the fit on the full data set we observe that the expansion up to the 
second order does not fit adequately the data: we find values of x^ probability very 
low, for any choice of the binning. Therefore we added to the amplitude expansion 
the cubic terms: 

\A{X,Y)\^ ^l + aY + hY'^ + cX + dX'^ + eXY + fY^ + gX^ + hX'^Y + lXY'^. (14) 

The values of probability improve significantly with the inclusion of only one 
additional term, namely the fY^ term; all the other cubic terms are consistent with 
zero and do not improve the fit quality; therefore in the following we keep only the 
/ parameter. 

Changing the bin size from AX = AY = 0.20 to AX = AY = 0.11 , we find very 
similar results, showing that the fit is almost binning independent. 
In table 1 are reported the fit results for AX = AY = 0.125 (154 bins are used), 
where we found the best x^ probability, and for different paramctcrisations of |Ap.In 
the first row all the parameters a, b, c, d, e, f are left free; in the following rows the 
C- violating c, e parameters are always set =0; in the last 3 rows we have moreover 



10 



set = 0; / = and = / = 0. 



dof 




10 3 a 


10 3 b 


10 3 c 


10 3 d 


10 3 e 


10 3 / 


147 


73% 


-1090±5 


124±6 


2±3 


57±6 


-6±7 


140±10 


149 


74% 


-1090±5 


124±6 




57±6 




140±10 


150 


< 10-6 


-1069±5 


104±5 








130±10 


150 


< 10-8 


-1041±3 


145±6 




50±6 






151 


< 10-6 


-1026±3 


125±6 











Table 1 

Fit results for different parameterisations of Values in 2"'^ row are used for final 
results. 

As expected from C-invariance in r/ ^ ti'^ti'tt^ decay the parameters c and e are 
consistent with zero, and removing them from the fit does not affect the result for 
remaining parameters. 

We observe a quadratic slope in X and a cubic slope in Y clearly different from zero. 
Therefore our final results for the Dalitz plot parameters are those shown in second 
row of the table. Fig. 5 shows a comparison between fit and data for the projections 
in X or Y. 




X Y 

Fig. 5. Comparison between data(points) and fit (histogram) for X,Y projections of the 
Dalitz plot distribution. 
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5 Systematic uncertainties 

Wc have estimated the systematics errors due to the following sources: 
efRciency evaluation All relevant reconstruction efficiencies have been checked 
directly on data, using suitable control samples: the only observable systematic 
effect is introduced hj the EVCL procedure, which we have check;ed with the 
minimum bias sample. We find a remark;able agreement between Data and MC 
for various kinematical distributions (see fig. 6); this implies that the experimental 
resolution is well modelled in MC; therefore we neglect any related systematics 
on the signal selection efficiency e {X, Y). 
resolution and binning We have checked energy resolution for the photons by 
comparing the distributions of the photon energies after the kinematic fit on data 
and MC finding good agreement for both the core and the tails of the distributions. 
Drift Chamber momentum resolution and absolute scale is controlled run by run 
by checking the Ks mass reconstructed in Ks 7r"'"7r~ events. The effect of 
binning was estimated by changing the bin size by a factor of two in the range 



background contamination The main source of backgrounds to our analysis are: 



(2) — * uJir^ with u — ^ tc'^tt'^tt^ ; 

Although overall background contamination is small, its presence gives an observ- 
able systematic effect at our level of precision. We have changed the cut on M^^ in 
a wide range, corresponding to B/S varying from 0.7% to 0.2% , obtaining slightly 
different values for the parameters, 
stability with respect to data taking conditions We have divided our data sam- 
ple in 9 periods of about 50 pb^^each. We have verified that the results for each 
parameter are compatible to be constant all over the data taking and that the 
average values over the periods are consistent with the ones from the whole data 
sample fit. 

radiative corrections Radiated effects have been considered by generating 10'' 
7] 7r+7r^7r°7 decays, according to ref. [16]. We have verified that the bin by 
bin ratio of the Dalitz plot obtained for 1] tt+tt^tt'^ 7 decays with the one for 
1] 71^71^71^ decays can be fitted with a constant with x^/dof = 154/153 corre- 
sponding to a probability of 46%. 



AX = AY e [0.11,0.20]. 
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For each effect mentioned above the systematic error has been conservatively esti- 
mated as the maximum parameter variation respect to the reference one; the total 
systematic error is obtained by summing in quadrature the various contributions as 
shown in table 2. 



Source 


Aa 


Ab 


Ad 


A/ 


EVCL 


-0.017 


0.005 


-0.012 


0.01 


binning 


-0.008 +0.006 


-0.006 +0.006 


-0.007 +0.001 


-0.02 +0.02 


background 


-0.001 +0.006 


-0.008 +0.006 


-0.007 +0.007 


-0.01 


Total 


-0.019 +0.008 


± 0.010 


-0.016 +0.007 


± 0.02 



Table 2 

Summary of the systematic errors on the Dalitz plot parameters. 
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Fig. 6. Data vs MonteCarlo comparisons in log scale. From top left in clockwise direction: 
minimum Pt and Pz, cosO between pion tracks and for photons. 
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6 Asymmetries 

While the polynomial fit of the Dalitz plot density gives valuable information on 
the matrix element, some integrated asymmetries are very sensitive in assessing the 
possible contributions to C violation in amplitudes with fixed A/. In particular 
left-right asymmetry (which is of course strongly related to the c parameter in our 
fit) tests C violation with no specific A/ constraint; quadrants asymmetry tests C 
violation in A/ = 2 and sextants asymmetry (for a definition see ref. [17]) tests C 
violation in A J = 1. 

In the following we present results on asymmetries which use 4 times the statistics 
entering the PDG fits. For this measurement we use the following approach: we 
obtain from MC the efficiency, for each region of the Dalitz plot, as the ratio between 
reconstructed and generated events in the region. This definition takes into account 
the resolution effects as well. We cross check the result by evaluating the asymmetries 
on Monte Carlo, these turn out to be all compatible with zero. We then evaluate 
the asymmetry on data by folding in the expected MC efficiency, after subtracting 
the MC expected background. On MC, for a sample of 5.69 x 10^ events we get: 

El = (34.91 ± 0.02)% Er = (35.05 ± 0.02)% Alr = (-0.006 ± 0.06) x 10"^ 
£13 = (35.01 ± 0.02)% £24 = (34.95 ± 0.02)% Aq = (-0.008 ± 0.06) x 10"^ 
£i35 = (35.00 ± 0.02)% £246 = (34.96 ± 0.02)% Ag = (-0.05 ± 0.06) x lO'^ 

The "raw" asymmetries on data are found to be: 

Alr = (-9 ± 9) x 10"^ = (2 ± 9) X 10"^ As = (13 ± 9) x 10"^ 

Correcting for MC efficiencies we get: 

Alr = (9 ± lOstat.) X 10-^; Aq = (-5 ± lOstat.) x 10"^ A5 = (8 ± 10,^^ x 10"^ 

In assessing the systematics for the measured asymmetries we have considered all 
the effects taken into account for the Dalitz plot fit, namely: 

• Effect of background (by varying cuts); 

• Effect of EVCL (by using Minimum Bias sample); 

• Data-MC efficiency comparison (using the (f) 7r"'"7r~7r° control sample) 

In particular the tracking efficiency has been evaluated separately for the two charges, 
since in MC is evident a small but statistically significant difference in left and right 
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efficiencies, due to a slightly different tracking efficiency as a function of px for 
positive and negative pions. 

Since we require both tracks to be reconstructed the absolute value of the efficiency is 
not important for the asymmetry, but rather its dependence upon the pion momen- 
tum. The very good Data-MC agreement has been already demonstrated for both 
charges on the signal. We here use an independent control sample o{ 6 ^ tt+tt^tt^ 
events and, as described before (see section 3), check the agreement between Data 
and MC for the efficiencies as a function of momentum the two charges. Results 
are shown in fig. 7. The control sample agrees well with MC within errors, and the 
Data-MC ratio is well fitted by a constant. In order to assess the possible systematics 
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Fig. 7. The Data-MC ratio of tracking efficiency for 7r+ (up) and vr^ (down) as a function 
of the pion transverse momentum. The maximum deviation from a constant fit has been 
used for systematics evaluation (see text). 

connected with the tracking efficiencies wc have assumed a "worst case" approach: 
we have estimated the maximum positive or negative linear slopes compatible within 
one sigma with the fit of the distributions shown in fig. 7. Then we have assumed the 
two charges to behave with opposite slopes. This gives us two possibilities (tt^ with 
positive slope and vr" with negative slope and vice- versa). Wc have then reweightcd 
the events according to these two possibilities and used the maximum difference 
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observed in the asymmetries as an estimate of the possible systematic effects due 
to different Data-MC efficiencies. The systematics connected with the asymmetries 
are shown in table 3. 



Syst. Effect 


Left-Right 


Quadrant 


Sextant 


Background 
EVCL 
Efficiency 
Total 


(-0.2/ + 0.1) X 10-^ 

(-0.5) X 10-=^ 
(-1.3/ + 0.9) X 10-3 
(-1.4/ + 0.9) X 10-3 


(-0.2/ + 0.2) X 10-3 

(-0.3) X 10-3 
(-0.3/ + 0.2) X 10-3 
(-0.5/ + 0.3) X 10-3 


(+0.3) X 10-3 
(+0.7) X 10-3 
(-1.3) X 10-3 
(-1.3/ + 0.8) X 10-3 



Table 3 



Systematic errors on asymmetries. 

Therefore the final results for the asymmetries are: 

Alr = (0.09 ± 0.10(stat.)in4(syst.)) x lO'^ 

Aq = (-0.05±0.10(stat.)l[J:[!^(syst.)) x 10"^ 
As = (0.08±0.10(stat.)l°:?^(syst.)) x 10-^ 

7 Conclusions 

The results including the statistical uncertainties coming from the fit and the esti- 
mate of systematics are: 

a = -1.090 ± 0.005(stoi)l°:°°^(sysi) (15) 

b = 0.124 ± Q.QQ6{stat) ± Q.QlQ{syst) (16) 
d = 0.057 ± 0.006{stat)1:oml(syst) (17) 
/ = 0.14 ± Om{stat) ± 0.02{syst) (18) 

The systematic error has been obtained adding in quadrature all the contributions 
in tab. 2. Tab. 4 gives the correlation coefficients between the fitted parameters. 
The following comments can be done: 
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a 


b 


d 


/ 


a 


1 


-0.226 


-0.405 


-0.795 


b 




1 


0.358 


0.261 


d 






1 


0.113 


f 








1 



Table 4 

Error matrix from the Dalitz plot fit. 

• the fitted value for the quadratic slope in Y is almost one half of the simple 
Current Algebra prediction {b — /^), thus caUing for significant higher order 
corrections; 

• the quadratic term in X is unambiguously different from zero; 

• similarly for the large cubic term in Y; 

• when integrating the polynomial over the phase space to get the decay width, the 
strong correlations between parameters must be carefully taken into account for 
a correct error estimate. 

• we don't observe any evidence for C violation in the 77 — > tt'^tt^tt^ decay since 
the c and e parameters of the Dalitz plot and the charge asymmetries are all well 
consistent with zero. 
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